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ABSTRACT 

Since its launch in late 2004, the Swift satellite triggered or observed an average of one gamma-ray burst (GRB) 
every 3 days, for a total of 771 GRBs by 2012 January. Here, we report the development of a pipeline that semi- 
automatically performs the data-reduction and data-analysis processes for the three instruments on board Swift 
(BAT, XRT, UVOT). The pipeline is written in Perl, and it uses only HEAsoft tools and can be used to perform the 
analysis of a majority of the point-like objects (e.g., GRBs, active galactic nuclei, pulsars) observed by Swift. We 
run the pipeline on the GRBs, and we present a database containing the screened data, the output products, and the 
results of our ongoing analysis. Furthermore, we created a catalog summarizing some GRB information, collected 
either by running the pipeline or from the literature. The Perl script, the database, and the catalog are available for 
downloading and querying at the HEASARC Web site. 
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1. INTRODUCTION 

Since the first gamma-ray burst (GRB) was observed by the 
Vela satellite (Klebesadel et al. 1973), the number of flashes 
of high-energy radiation that has been detected has increased 
dramatically, especially since the 1991 launch of the Compton 
Gamma Ray Observatory and its y-ray instrument, the Burst and 
Transient Source Experiment instrument (Meegan et al. 1992). 
Discoveries made by the BeppoSAX satellite starting in 1997 
showed that GRB afterglows are strong sources across a broad 
wavelength band (Sahu et al. 1997; Wijers et al. 1997). 

In 2004 November, NASA launched the Swift satellite 
(Gehrels et al. 2004), which has the ability to detect a GRB 
with its y-ray instrument (the Burst Alert Telescope, BAT; 
Barthelmy et al. 2005), to rapidly slew to the y-ray position, and 
to observe the GRB afterglow with the X-Ray Telescope (XRT; 
Burrows et al. 2005) and the Ultraviolet/Optical Telescope 
(UVOT; Roming et al. 2005). Through the end of 2012 
January, Swift observed 771 GRBs, an average of three GRBs 
detected every 10 days. The analysis of the data is performed 
by the Swift instrument teams, with each team focusing on their 
instrument. The results of these analyses are reported in the 
Gamma-ray Coordinates Network (GCN) circulars, a very ef- 
fective way of informing the GRB research community of the 
initial properties of GRBs. Although the circulars have an im- 
portant role, they also show some limitations. (1) The Swift 
data analysis usually refers to an earlier part of the GRB ob- 
servation (typically coinciding with the first few orbits), and 
only in a limited number of cases are the results from follow- 
up observations released. Furthermore, analysis methods have 
evolved and this may make some circulars partially incomplete, 
the results not uniform and, consequently, more difficult to com- 
pare. (2) The analysis in the GCN circulars is done with the 
software and calibration available at the time of the burst, and 
may not be updated when newer versions are released. These 
limitations suggested to us that an archive containing the re- 


sults of the data analysis of every observed GRB was desier- 
able. The following Swift instrument teams have processed data 
from their instruments: the NASA Goddard (GSFC) BAT Team; 
the UK Swift Science Data Centre at the University of Leices- 
ter for the XRT; and the UVOT team at Penn State University, 
NASA/GSFC, and the Mullard Space Science Laboratory. Each 
group focused on one specific instrument and a summary of 
their analysis is collected at the Swift Burst Ground-Analysis 
Web site (http://gcn.gsfc.nasa.gov/swift_gnd_ana.html), while a 
more complete database is available at the University of Leices- 
ter Web site (http://www.swift.ac.uk/xrt_live_cat/) for the XRT 
and at the “First Swift UVOT GRB Afterglow Catalog” Web site 
at NASA/GSFC (http://swift.gsfc.nasa.gov/docs/swift/results/ 
uvot_grbcat/) for the UVOT. The methods adopted by the BAT, 
XRT, and UVOT teams, as well as some of their results, have 
been extensively described (Sakamoto et al. 2008, 201 1 for the 
BAT; Evans et al. 2007, 2009, 2010, and Goad et al. 2007 for 
the XRT; Poole et al. 2008 and Roming et al. 2005, 2009 for the 
UVOT). 

In this paper, we present a pipeline that semi-automatically 
performs the data reduction and analysis of all the Swift 
instruments. We run the pipeline to analyze all the GRBs 
observed by Swift during its first 7 years of operation, and we 
describe the key products available for downloading. We also 
present our Swift GRB Catalog, whose purpose is to collect 
some basic information related to each GRB, retrieved from the 
analysis performed by the pipeline, the GCN Circulars Archive, 
and/or published papers. 

2. CATALOG: GCN INFORMATION RETRIEVAL 

After a GRB has triggered Swift or other X-ray/y-ray satel- 
lites like INTEGRAL , Fermi, AGILE, or the third Interplane- 
tary Network ( IPN ), the location of the burst is distributed by 
the GCN system. Within hours/days, the reports of follow-up 
observations made by ground-based and space-based optical. 
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infrared, radio, and X-ray observers are released through the 
circulars. If at any time the GRB is observed by Swift , the time 
of the detection, the most current and most precise position of the 
afterglow as well as possible information related to its distance 
(redshift) and the presence of a detected host galaxy are in- 
serted into our Swift Gamma Ray Burst Catalog (see Section 7), 
available at http://heasarc.gsfc.nasa.gov/cgi-bin/W3Browse/ 
w3catindex.pl. The trigger time and the position are the essential 
ingredients to run the Swift data analysis pipeline correctly. 

3. PIPELINE: DATA AND METADATA RETRIEVAL 

The first step in our automated pipeline consists of collecting 
all the files and metadata useful for the analysis of a specific 
source (e.g., GRB). Each observation is uniquely identified by 
an 1 1 digit number sequence: The first eight digits (the trigger 
ID) are associated with the GRB and the last three digits refer 
to different observations of that GRB. Using the trigger ID, the 
pipeline queries our Swift GRB Catalog for the time and po- 
sition of the burst, which are used as input for the script (see 
Section 4). The processing machine has direct, read-only access 
to the entire Swift archive, which allows the processing script 
to easily locate the sequence directories for each GRB. In each 
of these directories the processing script makes use of the un- 
screened event level 1 data for all the instruments used for the 
observations of a GRB as well as the UVOT level 1 raw images, 
if present. In addition, the pipeline uses the attitude and orbit 
data, the make-filter file, and the spacecraft attitude file from the 
auxiliary directory and the XRT housekeeping header packets, 
if the XRT was used for that observation. Only data obtained in 
the UVOT and XRT pointing modes are included in the analy- 
sis. For the XRT, these modes only include photo counting (PC), 
and window timing (WT) — piled-up mode data are not used (see 
Hill et al. 2004 for a detailed description of the XRT modes). 
PC mode is the typical frame transfer operation of a CCD and is 
used when the sources in the field of view (FOV) have low flux 
(since the readout time is 2.5 s, sources with a count rate higher 
than 0.4 counts s -1 will suffer a pile-up effect). Furthermore, 
PC mode retains full imaging and spectroscopic resolution. In 
contrast, WT mode loses some of the spatial information: this 
mode preserves only one-dimensional imaging by binning rows 
of the CCD and then reading out the central 200 columns of the 
CCD (corresponding to 8')- This enables high resolution timing 
(the readout time is 2.2 ms) with full spectral information for 
sources with very high count rates. 

For every other type of source, such as active galactic nuclei 
(AGNs), pulsars, etc., the user can specify the location of the 
directories that contain the level 1 data and the additional files 
required to run the analysis. 

4. PIPELINE: DATA REDUCTION 

The data reduction step of the pipeline is handled by a sin- 
gle script, swifthlprod (hereafter “the processing script”). 
This script performs both the data reduction and much of 
the analysis. It uses the latest release of the HEAsoft soft- 
ware and calibration database, available at the time of analysis. 
This processing script can be used with a standard 
HEAsoft distribution and is made available for download 
at http://heasarc.gsfc.nasa.gov/docs/swift/archive/grbsummary/ 
software.html. 

Before the processing script begins the data reduction in 
earnest, the attitude file for every observation is corrected for 
“jumps” in attitude data using the Swift task att jumpcorr. The 


jumps occur when the attitude is first updated using star tracker 
information after a slew. The corrected file is then used in the 
reduction of the XRT and UVOT data. 

For all the instruments, the processing script uses standard 
tools to produce the level 2 data: batgrbproduct for the BAT; 
xrtpipeline for the XRT; and coordinator, uvotmodmap, 
uvotimage, uvotbadpix, and uvotexpmap for the UVOT 
(currently there is not one single tool released by the instrument 
team that performs the complete data reduction). In the next 
sections we report in detail the set of parameters used for 
each tool by swifthlprod and how the user can control the 
output. 

4.1. BAT 

The standard BAT GRB processing script batgrbproduct 
is called by swifthlprod to do the data reduction of 
the event data, assuming the default values, i.e., aperture= 
“CALDB:FLUX,” short fix=“scaledmap, expand,” and a partial 
coding threshold of 5% (see http://heasarc.nasa.gov/lheasoft/ 
ftools/headas/batgrbproduct.html for a detailed description of 
these parameters). The processing script checks the value of the 
partial coding fraction reported in the output file report . txt, 
where a summary of the basic analysis results and infor- 
mation about the trigger are given. If this value is smaller 
than 0.05 (a threshold value indicating that the GRB posi- 
tion is poorly coded) or if batgrbproduct has failed to pro- 
duce the output directory and files, then the processing script 
runs batgrbproduct using the aperture parameter fixed at 
“CALDB DETECTION” and a minimum partial coding frac- 
tion of 0.0. If this trial fails as well, then batgrbproduct has 
to be run manually, assuming a case by case set of parameters. 
One example is when batgrbproduct fails to find a burst in- 
terval, as may be the case for short, weak GRBs. In this case, 
the operator must identify the trigger time and the stop time in 
order to create the standard level 2 products. 

4.2. XRT 

Nearly all the data reduction for the XRT PC and WT mode 
data consists of running the standard script xrtpipeline for 
each mode with the parameters set to the default values (e.g., 
for the PC mode, the threshold for the Pulse Height Analyzer, 
or PHA, of the central pixel is set to 80 and the minimum 
and maximum PHA values used are 0 and 4095, respectively). 
Only events with default grades of 0-12 for PC mode data 
and 0-2 for WT mode data are included in the cleaned event 
files. The processing script runs only the first two stages of 
xrtpipeline, i.e., it does not generate the standard scientific 
products. If xrtpipeline fails to produce an output cleaned 
event file for a particular mode, the data for that sequence / mode 
combination are not used. 

Occasionally, for PC mode data, xrtpipeline is unable to 
generate a uniform background across the CCD. This is typically 
due to either contamination by the South Atlantic Anomaly or 
Earth glow. To check for this, the following steps are taken: 
(1) the FOV is sub-divided into four quadrants, rotated to align 
with the PC chip, (2) the background level is estimated in each 
quadrant, and (3) the background in each quadrant is compared 
to the mean background of all quadrants. If the ratio of the 
mean background in any one quadrant to the mean background 
is larger than a specified threshold (currently set at 1.5), the 
pipeline will stall this sequence and await operator intervention 
(see also Puccetti et al. 2011 for another method to check for 
Earth glow). 
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4.3. UVOT 

Since there is no single UVOT data reduction tool provided in 
HEAsoft, swifthlprod performs the data reduction steps out- 
lined in the UVOT User’s Guide v2.2. The UVOT event mode 
data are minimally processed beyond the standard archive pro- 
cessing. The coordinator tool is run to update the SKY co- 
ordinates in the unfiltered event lists using the latest telescope 
definition and the corrected attitude files. The former contains 
information about the alignment of the UVOT instrument and 
the distortion of the image plane and resides in the calibra- 
tion database. The UVOT image data are checked for whether 
they have the modulo-8, correlated noise patterns removed. As 
explained in Section 2.1 of the UVOT User’s Guide and in 
Section 4 of Breeveld et al. (2010), the UVOT spatial resolution 
is increased because the algorithm on board Swift centroids the 
events generated by each incoming photon. Since the position 
of the centroid has a resolution eight times that of the original 
pixels, the software splits the pixels by a factor of eight to create 
a raw image. This process creates some artificial fixed-pattern 
structures that can be clearly seen in the raw images and have to 
be removed before performing accurate photometry. If any im- 
ages do not have this correction applied, then the uvotmodmap 
task is run on them. Both the image and event mode data are then 
combined into a single multi-extension FITS image for each 
filter/ sequence pair. In the case where an exposure has both 
event and image data, the former is not included in the analysis. 

These multi-extension FITS images are used as input for 
uvot image to produce sky images, and next used to calculate 
the aspect correction. The correction is estimated using, by 
default, at least five stars in the sky images. If the observation 
exposure time is too short or there is an insufficient number 
of bright stars, then the aspect correction is not calculated but 
must be interpolated from exposures that do have a good aspect 
correction. Only those exposures that are aspect-corrected are 
considered in the following analysis. The final steps of the 
UVOT data reduction are the running of the scripts uvotbadpix 
and uvotexpmap to create images containing the bad pixels 
and the exposure maps (for the latter, we leave the parameter 
method set to “MEANFOV”). As an additional output in the 
data preparation, the processing script adds all the exposures 
belonging to each observation in order to improve the signal-to- 
noise or to detect faint sources using the task uvotimsum. 

5. PIPELINE: DATA ANALYSIS 

All the steps described in this section are performed by 
swifthlprod. For each instrument, it calls individual tools that 
extract level 3 products such as light curves, spectra, and sky im- 
ages. The BAT analysis is done using a custom FTOOLS-like task, 
batgrbexprod. This task is also made available for download 
at http://heasarc.gsfc.nasa.gov/docs/swift/archive/grbsummary/ 
software.html. For the UVOT and XRT analysis, swifthlprod 
calls the HEAsoft tasks uvotgrblc and xrtgrblc. As ex- 
plained in Section 6, some of the level 3 products for the XRT 
require analysis tailored to the specific GRB. This is performed 
by the pipeline operator using the HEAsoft task xrtgrblcspec. 
Below, we give a detailed description of how these tasks perform 
the data analysis. 

5.1. BAT 
5.1.1. Light Curve 

Among the standard output of batgrbproduct is a set of 
light curves that cover the full range of event data in the entire 


BAT energy band 15-350 keV and in four sub-ranges (15-25, 
25-50, 50-150, and 150-350 keV). The time binnings of these 
light curves are fixed at 1 s, 64 ms, and 4 ms. In order to select 
the best time binning for that particular GRB, batgrbexprod 
estimates the length (L) of the observation (i.e., defined as the 
difference between the start and stop times of the total burst 
interval) and the highest intensity rate (r max ) observed in the 
64 ms light curve. Based on the combination of these two 
quantities, the best time binning (bin) is chosen as follows. 
(1) If r max ^ 1 then bin = 2 ms (L < 2.5 s), 64 ms (2.5 ^ L < 
100 s), and 256 ms (L ^ 100 s). (2) If r max < 1 then bin = 16 ms 
(L < 2.5 s), 256 ms (2.5 ^ L < 100 s), and 1 s (L ^ 100 s). 
Since batgrbproduct does not produce light curves with a 
binning of 16 and 256 ms, batbinevt is used to produce the 
light curves for the entire energy band and the four sub-ranges. 

5.1.2. Spectrum 

Spectra and response matrices are part of the standard output 
of batgrbproduct. Among others, spectra are computed for 
the total burst interval and the pre-slew, slew, and post-slew 
intervals. The response matrices are generated only for the pre- 
slew and post-slew intervals, since the BAT energy response 
generator, batdrmgen, works for a fixed single incident angle 
of the source, i.e., for the pre-slew and post-slew intervals, when 
the spacecraft is not slewing. For the slewing interval, we follow 
the method explained in Sakamoto et al. 2008: in summary, 
batgrbexprod (1) divides the interval into sub-intervals, each 
5 s long (if the last bin is shorter than 2 s, then it is appended to the 
second to last bin), (2) runs batbinevt to generate the PHA file 
corresponding to the bin, (3) updates the keywords of the PHA 
file (batupdatephakw), and (4) creates the response matrices 
with batdrmgen, taking into account the position of the GRB 
in detector coordinates. All the newly produced matrices are 
weighted and combined to create an average response matrix 
for the slewing interval using addrmf. The weighting factors 
for each matrix are estimated as follows: bin r /tot r , where bin, 
and tot c are the number of counts in each bin and in the entire 
slew interval. The spectra of the pre-slew, slew, and post-slew 
intervals are fitted within Xspec with two models, a simple 
power law and a cutoff power law (see Sakamoto et al. 2008 
for their definition). As explained in Sakamoto et al. (2008), a 
A/ 2 > 6 is the criterion to select the cutoff power law over the 
simple power law. 

A similar approach to creating files for a spectral analysis 
is adopted to convert the Bayesian block light curve into units 
of erg cm -2 s -1 . The conversion is obtained extrapolating the 
15-150 keV mask weighted (background subtracted) BAT light 
curve into the XRT 0.3-10 keV (and 2-10 keV) energy band. 
The selected light curve has Bayesian block bins and is part 
of the batgrbproduct standard output. This conversion com- 
pares the temporal source behavior covered by BAT and XRT 
observations over a long period of time. The conversion is based 
on the assumption that there is one continuous spectral behav- 
ior (a power law) over the 0.3-150 keV spectral range. Since 
the XRT and the BAT do not share a common energy range, 
this assumption had to be tested (Sakamoto et al. 2007). The 
choice of the 2-10 keV range is due to the fact that it is the 
XRT range closer to the BAT energy window. The Bayesian 
block input light curve can have bins that cross the time bound- 
aries of the pre-slew, slew, and post-slew intervals (e.g., part 
of the block is in one interval and the remaining in another 
interval). The PHA file and its response function correspond- 
ing to this bin has to be created accordingly, following the 
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Table 1 

XRT Extraction Regions 



PC 



WT 


Rate 

GRB 

BKG 

Rate 

GRB 

Pile-up 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

c/s > 7 

A(20",105") 

A(130", 260") 

c/s > 1000 

B(450") 

B(37'.'84) 

5 < c/s <7 

A(15", 95") 

A(120", 240") 

600 < c/s < 1000 

B(450") 

B(28'.'38) 

3 < c/s < 5 

A(12", 85") 

A(110", 220") 

400 < c/s < 600 

B(450") 

B(23'.'65) 

1 < c/s 3 

A(10", 75") 

A(100", 200") 

300 < c/s < 400 

B(300") 

B(9743) 

0.5 < c/s ^ 1 

A(5", 65") 

A(90", 190") 

100 < c/s < 300 

B(236") 

B(0") 

0.1 < c/s <0.5 

C(55") 

A(80", 185") 

50 < c/s < 100 

B(200") 

B(0") 

0.05 < c/s <0.1 

C(45") 

A(65", 165") 

10 < c/s < 50 

B(165") 

B(0") 

0.01 < c/s < 0.05 

C(25") 

A(40", 140") 

5 < c/s < 10 

B( 1 1 8") 

B(0") 

0.005 < c/s < 0.01 

C(22") 

A(40", 140") 

1 < c/s 5 

B(71") 

B(0") 

0.001 < c/s < 0.005 

C(17") 

A(40", 140") 

c/s 1 

B(35") 

B(0") 

c/s <0.001 

C(12") 

A(40", 140") 





Notes. Columns explanation: 1 = PC count rate range; 2 and 3 = PC source and background extraction regions; 4 = WT count rate range; 5 and 6 = WT source and 
pile-up extraction regions. The region shapes are: A = annulus, B = box, C = circle. 


method described above. For each bin, the spectral analysis is 
performed within Xspec, assuming a simple power law model. 
The count rate and fluxes in the 0.3-10 keV and 2-10 keV 
energy bands are estimated and used for the conversion. If a 
good fit cannot be found for a given Bayesian block, then this 
block is merged with the previous block, and the fitting pro- 
cess is repeated. The errors in each BAT light curve bin were 
found by randomly sampling the fitted model parameters within 
their respective confidence intervals 1000 times, and recalculat- 
ing the flux each time. The la confidence interval in this flux 
conversion sample is used as the error. Evans et al. (2010) de- 
scribe a somewhat different approach for combined XRT /BAT 
analysis, including a discussion of the many caveats that need 
to be considered when extrapolating from one instrument to 
another with a different bandpass. The results of this analysis 
for the growing list of Swift GRBs are available at the Web site 
http://www.swift.ac.uk/burst_analyser/. 

5.2. XRT 

Once the level 2 products are created following the processing 
described in Section 4.2, the task xrtgrblc determines the start 
and stop times of each orbit using the make-filter file stored in 
the auxiliary directory. This file contains selected spacecraft 
information related to the observation. The task considers as an 
orbit the time interval in which the XRT CCD is active (i.e., 
when the reference voltage for the variable “amp 1” is read 
from the filter file), the angular distance from the target is less 
than 10', the spacecraft is settled on the target, and when the 
gap between two consecutive X-ray photons reaching the CCD 
is less than 100 s. At this stage the data preparation follows 
one of the following two paths, depending on the mode of the 
observation. 

1. WT mode. Despite having selected only times where the 
spacecraft is considered settled, some small adjustments in 
the pointing can happen, in particular, at the very beginning 
of the observation, when the observation mode is most 
likely to be WT. The orbit boundaries are modified in order 
to consider only those times when the spacecraft roll angle 
is less than 1° from the final settle point. 

2. PC mode. Before the light curve extraction and analysis 
begins, the location and intensity of any field sources in the 
FOV are determined using Ximage. These two quantities 


define the size and location of regions that must be excluded 
from the light curve analysis. 

5.2.1. Light Curve 

To have a fully corrected XRT light curve, adjustments have 
to be considered to address some issues, e.g., the extreme source 
variability, the pile-up effect for bright sources, and some CCD 
artifacts (Evans et al. 2009). The variability and the pile-up 
effects can be taken care of by selecting the correct size of 
the source, background, and pile-up extraction regions. The 
selection is dependent on the event data mode, and the two 
modes (PC and WT) each warrant their own discussion. For 
the PC mode a preliminary light curve for each orbit from a 
40" circle centered on the source is extracted using the FT00LS 
extractor and with a time binning of 20 s. A sliding-cell-type 
algorithm is applied to this preliminary light curve to determine 
the time-dependent source and background extraction regions, 
and any pile-up exclusion. This algorithm uses the running 
weighted mean count rate to divide the light curve into time 
intervals during which the count rate falls within empirically 
determined ranges. Each count rate range has a corresponding 
source and background region shape. Using this procedure, 
during periods of high count rate, larger source and background 
regions are selected. The regions for the source are either circular 
or annular (the latter to avoid pixels affected by pile-up), while 
the background regions are always source-centered annuli. The 
region size for the source (GRB) and background (BKG) are 
reported in the left section of Table 1 . 

The source regions were selected to maximize the contribu- 
tion of the source and to minimize the number of background 
photons. A reverse approach has been adopted to determine the 
background region size, with the additional caveat of consider- 
ing areas relatively close to the source. Circular regions associ- 
ated with field sources are subtracted from both the source and 
background regions. Their radii correspond to twice the half- 
box length reported by Ximage. If the field source is within 0'6 
of the source, then this radius is reduced by half. 

Sometimes the source position is close to the edges instead 
of the center of the CCD. For GRBs, this will typically happen 
when the detection position is not well known, e.g., when the 
GRB is discovered by missions that carry instruments with poor 
spatial resolution, like Fermi or the IPN. In this case, the size of 
the source extraction region will be reduced from the empirical 
values of Table 1 in order to be fully contained in the CCD. 
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The background region will then be a circle of 150" located 
close to the source, fully contained in the CCD, and with field 
sources excluded. 

For the WT mode event data, both the preliminary and 
final light curve extraction regions are positioned over the WT 
strip using a combination of the HEAsoft tools pointxform, 
sky2xy and xy2sky, and the spacecraft housekeeping data. The 
regions themselves are all 20" wide. Similar to the PC mode, 
a preliminary light curve is extracted for each orbit from a 50" 
long box centered on the source using the FTOOLSextractor 
and with a time binning of 4 s. The above sliding-cell algorithm 
is applied to this preliminary light curve. As shown in Romano 
et al. (2006), the pile-up effect can be present for very bright 
sources. Following their results, we selected as source and pile- 
up regions, boxes centered on the source position and with 
length as a function of the counts rate. The sizes of these boxes 
are reported in Table 1 . The background regions are a 450" long 
box centered on the WT strip, excluding the source region plus 
an additional 55" on either side to avoid point-spread function 
(PSF) contamination. Field sources are ignored in the WT mode. 

The task xrtgrblc extracts light curves in the 0.3-10 keV, 
0.3-2 keV, and 2-10 keV energy bands. These curves are 
uniformly binned and the bin duration used is 0.01766 s for the 
WT mode and 2.51 s for the PC mode. The source light curve 
must be corrected for the presence of pile-up, irregularities in the 
exposure map (e.g., bad CCD pixels and columns, and attitude 
variations), and the PSF. The correction factor is estimated using 
the HEAsoft tool xrtlccorr. The background light curve is 
corrected only for irregularities in the exposure map using the 
HEAsoft tools xrtexpomap and Ximage. The output of this is 
subtracted from the raw source light curve and then multiplied 
by the total source correction factor to give the net corrected light 
curve. The total PC and WT corrected light curves are binned 
using a step-binning algorithm, where the number of counts 
included in each bin is proportional to the average rate in that 
bin. The algorithm requires that each bin has at least 20 counts 
and a signal-to-noise ratio of at least 3. Often at the end of the 
PC light curve, there are bins that do not satisfy these conditions. 
In this case, a 3a Bayesian upper limit is calculated. Since the 
WT mode is active only for high-count-rate epochs, low signal- 
to-noise bins are not converted to upper limits and are discarded 
from the binned light curve. This method is also applied to the 
hard (2-10 keV) band. The hard band bin boundaries are applied 
to the soft (0.3-10 keV) band, and a hardness ratio, defined as 
hard vs. soft, is calculated. 

A source and background spectrum for each of the time- 
dependent extraction regions found previously are extracted. 
The background-subtracted source spectrum in the 0.3-10 keV 
energy range is fitted with a doubly absorbed power law. One 
of the two absorption components is fixed to the Galactic value 
and the other is left free to vary. This spectral fit is necessary to 
obtain the flux conversion factors in the 0.3-10 keV, 0.3-2 keV, 
and 2-10 keV energy ranges. If the fit is good (0.7 ^ / 2 ^ 1.3), 
the conversion factor is estimated as the ratio of the model flux 
to the count rate. If not, then the flux conversion factor from the 
previous interval is used. If there is no prior good fit, then mean 
source flux conversion factors of 5 x 10 _n erg cm -2 count -1 
(0.3-10 keV), 2.7 x 10 -11 erg cm -2 count -1 (0.3-2 keV), and 
2.3 x 10 -11 erg cm -2 count -1 (2-10 keV) are used. 

The script xrtgrblcspec calls Xspec to fit the corrected 
0.3-10 keV light curve with three models that reproduce the 
typical trends observed in GRB light curve decay: a simple, a 
broken, and a double broken power law. The best fit is found 


using a series of f -tests, with a significance threshold of 3a . If 
there is some GRB flaring activity, the above models will not 
adequately represent the light curve. For this reason, we always 
visually inspect the light curve fits and manually remove flares 
if necessary. By doing so, we ensure that the fit only reproduces 
the continuum light curve. The human intervention also ensures 
that the fit is as optimal as possible. There are a few cases 
where the light curve should be fitted with a more complex 
model than the double broken power law; Evans et al. (2009) 
found that 110 light curves out of 665 required more than two 
breaks of 4a. Unfortunately, Xspec does not provide a power 
law model with a variable number of breaks. In those cases, the 
user tries to minimize the residual by using the double broken 
power law model. 

5 . 2 . 2 . Spectrum 

The breaks in the XRT light curve can emphasize a change in 
the environment where the X-rays are emitted or a change in the 
physical properties of the emitting region. We use these breaks 
to define time intervals where we expect to find no spectral 
variability. If the light curve is best fitted with a simple power 
law, the entire observation is considered as a single interval, 
xrtgrblc automatically sums all the source 0.3-10 keV PHA 
files and weighted ancillary response functions (ARFs) for the 
orbits that belong to each time interval in the light curve. The 
ARFs are generated using the xrtmkarf task (Moretti et al. 
2005). The sum is obtained running the FT00LS tasks mathpha 
and addarf (the error method used when calling the mathpha 
task is Gaussian). The reason for summing the PHA file is 
twofold: (1) the number of orbits can be large enough to create 
problems for Xspec in fitting each orbit as a single data set 
and (2) the better signal-to-noise makes it easier to see spectral 
features in the combined spectrum than in the individual spectra. 
The combined spectra are fitted using a double absorbed power 
law, where one absorption is fixed at the Galactic value and 
the other is left free to vary. The spectral analysis can reveal 
a change in the level of absorption or in the intrinsic photon 
distribution. 

5.3. UVOT 

The task uvotgrblc uses uvotsource to extract photomet- 
ric light curves from UVOT images. The accuracy of these 
curves is ensured by applying the coincidence-loss and large- 
scale sensitivity correction (estimated by uvotsource) and by 
accounting for several known issues: (1) presence of contami- 
nating field sources, (2) selection of the background extraction 
region, (3) significance of the detection, (4) position, and (5) size 
of the source aperture. These issues are addressed as follows. 

1. Field sources and background region. The choice of the 
background region depends on two factors: the presence 
of field sources (i.e., their intensity and their location with 
respect to the source) and the intensity of the source (i.e., 
how far the PSF wings extend). The first step is to run 
uvotdetect to obtain a list of field sources and their 
exclusion regions in each UVOT summed image. Next, 
uvotgrblc selects the optimal region size by estimating 
the background of the summed images in three annular 
regions centered on the position of the source (in the case 
of our analysis of GRBs, it is derived from our Swift GRB 
Catalog). The annuli have inner-outer radii of 15-27", 
27-35" , and 35^-2", respectively, and detected field sources 
are excluded. We are confident that the choice of these 
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radii is conservative and no significant contribution of the 
source PSF is expected, at least in the largest annulus 
(see Figure 5 in Poole et al. 2008, where it is shown that 
the wings of the PSF can be seen up to a maximum of 
25" from the source center). The three background values 
are compared to each other and the region with the lowest 
value is selected. There is a chance that all three regions 
were contaminated by photons from a nearby source (e.g., 
a galaxy) or CCD artifacts (e.g., ghost images or readout 
streaks). The background region is visually inspected for 
these types of contaminants. If there is a problem, a custom 
background region is selected and the processing script is 
run again for this source. 

2. Significance of detection. An approximate significance is 
initially determined by extracting the number of counts in 
a 3" radius circle and in the background region (see above) 
for all the exposures in every filter. For dim sources, the 
source aperture can be background-dominated. Thus, we 
consider as a detection every time Odet > 3. We use fT det = 
counts ne ,/(T to t, where we define cr tot = (a s 2 0urce + a b 2 kg ) 1/2 . In 
the estimate of counts net , we use the total corrected counts 
if the counts frame” 1 is >1, and the total corrected counts 
scaled by a factor of 

total counts 
number of frames 

otherwise. This is a conservative approach that does not 
allow us to miss weak sources. 

3. Aperture position. The extraction region must be centered 
exactly on the source. If there is a detection in any filter 
(see above), the coordinates of the object closest to the 
position provided by the user (or by the Swift GRB Catalog 
in our analysis) are refined using the task uvotcentroid. 
The task is run on the exposures with the highest detection 
level for each filter. The average right ascension (R.A.) and 
declination (deck) are found after a single round of la 
rejection of outliers. If the new position is beyond 3" from 
the catalog position, then no adjustment is made and the 
catalog’s R.A. and deck are used. If the new coordinates 
are offset from the inputs by more than 1"5, then the initial 
estimates of the significance of detection in each exposure 
are recalculated, as well as the optimal background regions. 

4. Aperture size. The standard photometric aperture is a circle 
with a 5" radius (Poole et al. 2008). This area is large enough 
to collect most photons for a moderately bright object. For 
weaker objects, the PSF wings are too dim and the area 
is mostly background-dominated: the radius of the source 
extraction region must change accordingly to determine the 
significance of detection and to extract the best photometric 
information. The intensity and position of the field sources, 
combined with the initial source intensities, are used to 
determine the optimal size for the photometry aperture. 
The selection of the aperture radius is as follows: (1) 5" if 
the nearest neighbor is distant (>15") or is closer (>10") 
but dim (<0.1 c/s), (2) 4" if the nearest neighbor is closer 
(>10") but brighter (0.1 < c/s ^ 10), (3) fixed at 3" for 
all other cases. The radius can be reduced from 5" and 
4" to smaller dimensions based on the source intensity, if 
detected: for c/s <0.1, the radius becomes 3"; for 0.1 < 
c/s < 0.75 the radius is fixed at 4". 

If the source extraction region is smaller than 5" then the 
photometric estimates derived using uvotsource have to 
be aperture-corrected. To estimate the correction factor, the 


three “best” field sources for every OBS_ID are selected. 
Each source must meet the following selection criteria to 
be considered a good candidate for aperture correction: (1) 
be nearly circular in shape according to uvotdetect (i.e., 
the ratio of the semimajor to semiminor axes is less than 
1.5), (2) lie in the range of 20" to 200" from the source, 
(3) have a size <7", (4) have a count rate in the range of 
0.5 c/s to 5 c/s, (5) not have any nearby neighbors, and 
(6) have an average exposure >95% of the summed image 
exposure. The ratio of the count rates obtained in circular 
regions with 5" radius and with 3" or 4" gives the aperture 
correction factor for each of the three stars. The pipeline 
selects as optimal the median of the three values, limiting 
the possibility of error in the selection of the three objects. 
If three good sources cannot be found by this method, then 
the aperture correction is estimated using the parameter 
apermeth=CURVE0FGR0WTH in the task uvotsource. 

If the source was not detected in any consecutive images, 
those images are summed and photometry is performed on the 
result. If in any image, or summed image, no source is de- 
tected, an upper limit is estimated. The photometric measure- 
ments and the upper limits are derived from the FLUX_HZ 
(and FLUX_HZ_ERR) and FLUX_HZ_LIM columns, respec- 
tively, in the uvotsource output FITS file (obtained with the 
parameter apercorr set to NONE). They are expressed in units 
of flux density and are not suitable to generate a combined light 
curve in which all the filters are compared. All these quantities 
can be converted to units of flux (erg cm” 2 s” 1 ). The formulas 
used to estimate the flux for the detection, its error, and for the 
upper limit are the following: 

FLUX D et = FLUX_HZ x AREA x APERTURE x IE-26 
FLUX_ERR D et = FLUX_HZ_ERR / FLUX_HZ x FLUX Det 

FLUX Up = FLUX_HZ_LIM x AREA x IE-26, 

where APERTURE is the correction factor, and AREA is the 
integrated filter response convolved with a power-law spectrum 
with index a = 1 . The functions that reproduce the filter response 
can be found in the CALDB files swuXX_20041 120vl0Y.arf, 
where XX indicates the filter and Y the version. Since all the 
flux densities are multiplied by the filter response, the integrated 
flux light curve will present the white light curve above all the 
other filters, due to its larger band pass (e.g., see Figure 12, right 
panel). 

6. PIPELINE: OUTPUT PRODUCTS 

From the level 3 products (light curves, spectra, images) and 
the results of the temporal and spectral analyses, a series of Web 
pages are automatically generated for each source. The typical 
products for the three instruments are gif images and FITS files. 
The images show a small area of the sky around the source direc- 
tion, the light curves, and the plots of the spectral analysis (for 
the BAT and the XRT only). The FITS files can be divided into 
three main categories: light curves, spectral files (PHA, ARF, 
and response matrix files or RMF), and a summary file for the 
BAT and XRT instruments. The summary file contains all the 
results of the light curve and/or the spectral analysis, together 
with the indication of important quantities used in the data anal- 
ysis. These quantities are used to populate various tables con- 
tained in the Web pages. The files produced by our pipeline are 
available for download from the HEASARC archive Web site: 
http://heasarc.gsfc.nasa.gov/docs/swift/archive/grbsummary/ 
grb_pages/20XX/grbXXXXXXX/. For an explanation of the 
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GRB09061 8: BAT 15-350keV, TARGJD 355083 


Pre— Slew, Exposure 59.0s 



Post— Slew, Exposure 852.8s 
22 h 20 m 00 s 21 h 40 m 00 s 20 h 20 m 00 s 18 h 00 m 00 s 



Figure 1 . BAT 15-350 keV FOV in the pre-slew and post-slew intervals. 
(A color version of this figure is available in the online journal.) 


+85 00 00 


+80 00 00 


+75 00 00 


GRB09061 8: BAT, TARGJD 355083, Bin time=0.256s, 1 Channel 



Figure 2. BAT 15-350 keV background subtracted light curve. 
(A color version of this figure is available in the online journal.) 


naming convention adopted for the files, see the document 
at http://heasarc.gsfc.nasa.gov/docs/swift/archive/grbsummary/ 
Filenames.pdf. A detailed description of all the available prod- 
ucts is given below. 

6.1. BAT 

6.1.1. Images 

The script batgrbproduct generates three sets of gif images 
for the BAT instrument and they provide visual information 
of the spatial, temporal, and spectral analyses. The three sets 
contain the following: 

1 . The FOV of the detected burst derived from the event data 
during the pre-slew and, if available, after the slew centered 
on the GRB position. Both images, standard products of the 
task batgrbproduct, show the source in the 15-350 keV 
energy range and are background-subtracted (Figure 1). 

2. The mask-weighted light curves of the full GRB duration 
(T100, that is, the time interval from the emission start 
time to the emission end time). T100 is determined by 


battblocks using a Bayesian-block technique (Sakamoto 
et al. 201 1). The light curves are obtained with the optimal 
time binning in the total 15-350 keV light curve (Figure 2) 
and in four energy ranges, corresponding to 15-25, 25-50, 
50-100, and 100-350 keV (Figure 3). Both plots have an 
additional 5% of the T100 interval at both boundaries for 
plot purposes. The y-axis unit is counts s -1 det~\ where 
a det is 0.4 x 0.4 = 0.16 cm 2 , and the time binning is 
chosen based on the length and intensity of the burst (see 
Section 3). In the second plot, the slew interval, if present, 
is marked with a horizontal line. 

3. The results of the spectral fit with a power law and cutoff 
power law in the 15-150 keV range for each epoch, i.e., 
pre-slew, slew, post-slew, and total (Figure 4). 

6.1.2. FITS 

The FITS files created by batgrbproduct, used to generate 
the images described above and to summarize the spectral 
analysis, are the following: 


7 




The Astrophysical Journal Supplement Series, 203:2 (17pp), 2012 November 


Donato et al. 


in 

CM 

I 

in 


o 

m 

I 

m 

CM 


o 

o 


I 

o 

m 


> 


Q) 

o 

in 

m 

I 

o 

o 



0 50 100 150 200 

Time Since BAT Trigger (s) 

Figure 3. BAT background-subtracted light curves in the four energy ranges (15-25, 25-50, 50-100, and 100-350 keV). 
(A color version of this figure is available in the online journal.) 



Figure 4. Spectral fit of the BAT pre-slew, slew, post-slew, and the total intervals with the simple power law (top) and power law with cutoff (bottom). 

(A color version of this figure is available in the online journal.) 


1. The fixed time bin standard one- and four-channel light 
curves as produced by batbinevt. 

2. The Bayesian binned light curves (see Section 5. 1) in units 
of erg cm~ 2 s” 1 in the ranges 0.3-10, 2-10, 15-25, 25-50, 
50-100, 100-150, and 1 5-1 50 keV. For each bin, the values 
of the spectral fit and the statistic are given. 

3. A PHA spectrum and corresponding response matrix for 
each epoch (i.e., pre-slew, slew, post-slew, and total). 

4. The information file (BATINFO extension), containing all 
the values of the spectral analysis, e.g., spectral index of 
the power law (with errors) and its normalization, cutoff 


energy (with errors), model flux in various energy ranges, 
degrees of freedom, and / 2 3 4 value. 

6.2. XRT 
6.2.1. Images 

Similarly to the BAT, there are three sets of gif images for the 
XRT. The script xrtgrblc generates the following: 

1. The 0.3-10 keV images centered on the source with a 
large and a small FOV. In both modes the plots show the 
extraction regions for the source (in solid blue) and the 
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30* 37 m 00* 30* 19 h 36 m 00* 30* 35 m 00* 34 m 30* 



GRB09061 8: XRT 0.3-10.0 keV 


TARGJD 355083 



Figure 5. Top: 0.3-10 keV WT mode large and small FOVs of the first observation. Bottom: 0.3-10 keV WT mode small FOV of all the orbits for all the observations 
combined. 

(A color version of this figure is available in the online journal.) 


background (in dashed blue) and the excluded areas (in 
dashed red), such as the field sources in the PC mode, and 
the pile-up region in the WT mode). The PC image is the 
sum of all the observations (Figure 5, upper panel), while 
the WT image refers to the interval with the highest signal- 
to-noise ratio (Figure 6, upper panel). 

2. A mosaic of the images of each single observation ID, to 
show the detection and variability of the source during the 
entire observational campaign (Figure 5 and Figure 6, lower 
panels, for WT and PC mode, respectively). 

3. The 0.3-10 keV combined WT and PC binned light 
curves in units of count rate (Figure 7, upper panel). 
The position in time when the shape of the source ex- 


traction regions change are indicated by up arrows. The 
blue and red colors are adopted for the WT mode and 
the PC mode, respectively. At the top of the plot the 
OBS_ID time boundaries are shown. The bottom part of 
the plot reproduces the fractional exposure of each light 
curve bin. 

4. The 0.3-10 keV combined WT and PC binned light curves 
in units of erg cm~ 2 s _I (Figure 7, lower panel). The 
procedure to estimate the conversion factors from rate to 
flux is explained in Section 5.2. 

5. The combined WT and PC hardness ratio, defined as the 
ratio of the 2-10 keV to 0.3-2 keV binned light curves 
(Figure 8). 
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Figure 6. Top: 0.3-10 keV PC mode large and small FOVs of the first observation. Bottom: 0.3-10 keV PC mode small FOV of all the orbits for all the observations. 
(A color version of this figure is available in the online journal.) 


6. The combined WT and PC 0.3-10 keV binned light curves 
in units of count rate with the best model obtained during 
the light curve analysis (Figure 9). The dashed black lines 
indicate the position at the time of the light curve breaks. 
Sometimes the fitted values for the time breaks may not 
exactly match the spectral intervals fitted on the spectra 
page. This is because the time boundaries used to derive the 
spectra correspond to the nearest light curve bin boundary. 
In this Xspec-style plot, the bottom part reproduces the 
residual of the data versus the model. 

7. The combined WT and PC 0.3-10 keV spectrum of each 
interval, defined during the light curve analysis (Figure 10). 
These spectra refer to the continuum emission of the source, 
as explained in Section 5.2. 

6.2.2. FITS 

The input FITS files for the light curve and spectral analyses 
and the summary of the results of both analyses can be 
downloaded. The user will find the following: 


1 . The fully corrected and background-subtracted PC and WT 
binned light curves, in the entire energy range of 0.3-10 keV 
and in the two sub-ranges of 0.3-2 and 2-10 keV. These 
curves are both in units of count rate and erg cirT 2 s -1 . 
Together with them, the source PC and WT hardness ratio 
is also present. All these curves have been used to generate 
the images previously described. 

2. The fully corrected PC and WT unbinned light curves of 
the source and the background in the entire energy range 
0.3-10 keV and in the two sub-ranges 0.3-2 and 2-10 keV; 

3. PH A spectra for the source and the background and 
the corresponding response matrix for one or more time 
intervals and for both the WT and PC modes. The time 
intervals were determined with the light curve analysis, as 
described in Section 5.2. 

4. The information file: XRTINFO, XRTFLDREG, and 
XRTF1TINFO extensions. XRTINFO contains the name 
of the files used for the data reduction and analysis, in- 
formation used to create the final light curve for each of 
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Figure 7. Combined PC (blue) and WT (red) mode binned light curve in unit of count rate (top panel) and of erg cm 2 s 1 (bottom panel) in the 0.3-10 keV energy 
range. The fractional exposure shown below the count rate light curve is defined as the exposure time divided by the bin duration. 

(A color version of this figure is available in the online journal.) 


the orbits, such as the start and stop times, the source 
and background extracted counts (in the 0.3-10 keV, 
0.3-2 keV, and 2-10 keV ranges), the spatial and tem- 
poral correction factors, the count rate to flux conver- 
sion factor, and details of the source and background ex- 
traction regions. XRTFLDREG collects all the data for 
the detected FOV sources found in the PC event file by 
Ximage: position, intensity, and probability of detection, 
as well as the size of the region associated with each 
source. This information is necessary to create the series 
of exclusion regions that are used to extract the source 
and background light curves and spectra (see Section 5.2). 
XRTFITINFO contains the results of both the light curve 


and spectral analyses. For every fit, it reports the times 
of the interval that has been fitted, the mode (WT or PC) 
used to collect the data, the fit statistics (degrees of free- 
dom and x 2 value), which model has been adopted (in 
Xspec style), and who performed the fit (in automatic by 
the pipeline or by the intervening user if the automatic fit 
failed or was not optimal). Values up to six parameters 
are reported. These parameters refer to different quanti- 
ties, based on the type of fit, e.g., light curve or spectral 
fit: (1) first slope or Galactic absorption, (2) first break 
or intrinsic absorption, (3) second slope or spectral in- 
dex, (4) second break or spectral model normalization, (5) 
third slope, (6) light curve model normalization. For the 
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GRB090618: TARGJD 355083, XRT Hardness Ratio 
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Figure 8. Hardness ratio for the PC (blue) and WT (red) data. The hardness ratio, defined as H/S, where S and H are the number of counts in the soft (0.3-2 keV) and 
hard (2-10 keV) bands, respectively. The count rate light curves in the two bands are plotted in the top and central panels, while the hardness ratio is plotted in the 
bottom panel. We adopted the binning of the hard band light curve to group the soft band light curve. 

(A color version of this figure is available in the online journal.) 


spectral analysis, the value of the observed flux is given. All 
these quantities come with the upper and lower value, de- 
termined after adding and subtracting the related ler error, 
respectively. 

6.3. UVOT 

Since no spectral analysis is performed on the UVOT data, 
only two sets of gif images are generated by uvotgrblc, 
providing spatial and temporal information. The sets include 
the following: 

6.3.1. Images 

1. The summed image with the highest signal-to-noise ratio 
centered on the source and with an FOV of ~4!5. For 
this particular OBS_ID and filter, the aperture used for 
photometry, the background region, and the excluded field 
sources are indicated with a solid blue circle, a blue dashed 
region, and red dashed circles, respectively (Figure 11, left 
panel). 

2. The earliest exposure centered on the source and with an 
FOV of ~7I5. Overplotted on this image are the BAT 
and XRT error circles (Figure 1 1, right panel). For GRBs, 
sometimes this exposure is obtained when the satellite is 
still settling. This will result in sources with elongated trails, 
as in the case presented in Figure 1 1 . 

3. The light curves in units of magnitude and erg cm -2 s -1 
(Figure 12). For each filter, the conversion factors from flux 


density to integrated flux have been estimated as explained 
in Section 6. In the plots, down arrows indicate the 3a 
upper limit for that bin. The colors in the light curve are in 
units of erg cm -2 s -1 to match those used in the magnitude 
light curve: red = V, green = B, blue = U, light blue = 
Wl, magenta = M2, yellow = W2, orange = WHITE. 
We want to point out that the light curves may contain the 
first exposure in the V filter, which can be affected by the 
trailing issue (see the previous bullet). This will reduce the 
real count rate and, consequently, will affect the estimate of 
the magnitude and the flux. We urge the user to discard this 
point in the light curves if the trailing effect is seen in the 
image, as in Figure 11, and to regenerate the light curves 
with a customized analysis. 

6.3.2. FITS 

The only UVOT FITS file available for download is the input 
file used to generate the two light curve plots (see above). 
For each exposure or sum of exposures it reports the final 
photometric value for the detection or the upper limit. It is 
produced by uvotsource, but is modified by uvotgrblc to 
contain additional information, such as the integrated flux (with 
error), the estimated sigma detection, the inner and outer radii 
of the background region, the aperture used to extract the 
source photometry, the aperture correction factor, and the basic 
information about the star used to determine this factor. As in 
the case of the light curve plots, the FITS file may also contain 
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GRB090618: XRT PC/WT Automatic Light Curve Fit 
Double Broken Power— Law 



Figure 9. XRT light curve analysis. The dashed lines indicate the position of the breaks found by the script xrtgrblcspec, using the double broken power law model 
within Xspec. The lower panel shows the ratio obtained by dividing the data by the folded model to highlight deviations. 

(A color version of this figure is available in the online journal.) 


the value from the first exposure in the V filter. This should be 
ignored if Swift was still settling at the time of the data collection 
the user must perform a customized analysis to obtain reliable 
results. 

6.4. Combined XRT and BAT 
6.4.1. Images 

The processing script swif thlprod also generates two plots 
that show the combined XRT and BAT light curves in the 0.3-10 
and 2-10 keV energy bands (Figure 13). The Bayesian block 
BAT light curve (in black) is combined with the WT (in blue) 
and PC mode (in red) binned light curve. The BAT light curve 
conversion to these two energy bands is explained in Section 5. 1 . 

7. CATALOG: PIPELINE AND ADDITIONAL 
INFORMATION RETRIEVAL 

Once the GRB analysis is completed, some quantities related 
to the XRT and UVOT analysis are extracted from the files 
generated by the pipeline and imported in the Swift Gamma Ray 
Burst Catalog (see Section 2). 

For the XRT, this table provides information regarding the 
following. 

1 . The first 100 count interval: the corrected 0.3-10 keV count 
rate, the start and stop time, exposure and mode used (PC, 
WT, or both). 


2. The first 300 s interval: the corrected count rates in the 
total 0.3-10 keV band, as well as the three sub-bands 
(0.3-1 keV, 1-2 keV, and 2-10 keV), the start and stop 
time, exposure, mode used, and the hardness ratios between 
the sub-bands. The rate is estimated as follows: (a) if there 
are at least 100 counts detected across 300 s, then the time 
range corresponds to the arrival of these counts, (b) if there 
are less than 100 counts, then the time range is the full 
observation (if there is not a detection, an upper limit is 
provided), (c) if the observation starts 400 s or more after 
the GRB trigger time, then no value is reported. 

In addition, the table provides the presence of flares and 
light curve changes in the continuum and the most accurate 
position observed with the XRT. The latter is derived from the 
UK Swift Science Data Centre Web site at the University of 
Leicester (http://www.swift.ac.uk/xrt_positions/). The position 
is typically enhanced by correcting the astrometric reference 
of the XRT. This will reduce the absolute systematic error 
on the position to 1"4. The correction is obtained by “matching 
the stars seen by UVOT with the USNO-B1 catalog and using 
the known mapping between the UVOT detector and the XRT 
detector.” 7 A detailed description of the method used can be 
found in Goad et al. (2007) and Evans et al. (2009). 


1 Swift Science Data Centre Website, http://www.swift.ac.uk/xrt_positions/ 
docs.php. 
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GRB09061 8: XRT PC/WT First Interval Spectra (0.3-10.0 keV) 
Interval [T 0 + s]: 124.7 - 306.3 

Blue - WT Mode 


GRB090618: XRT PC/WT Second Interval Spectra (0.3-10.0 keV) 
Interval [T 0 + s]: 306.3 - 19986.0 

Red - PC Mode Blue - WT Mode 



0.5 1 2 5 

Energy (keV) 

GRB090618: XRT PC/WT Third Interval Spectra (0.3-10.0 keV) 
Interval [T 0 + s]: 23295.4 - 2808500.0 

Red - PC Mode 


Energy (keV) 

GRB090618: XRT PC/WT Total Spectra (0.3-10.0 keV) 
Interval [T 0 + s]: 124.7 - 2808500.0 

Red - PC Mode Blue - WT Mode 




1 2 
Energy (keV) 


Figure 10. Plots of the XRT spectral analysis at different time intervals. The time breaks were obtained by fitting the light curve with a double broken power law 
model in Xspec (see Figure 10). As in previous figures, the blue and red points refer to WT and PC mode data, respectively. 

(A color version of this figure is available in the online journal.) 


For the UVOT, our Swift Gamma Ray Burst Catalog reports: 
(1) the observed magnitude and monochromatic flux of the first 
observation in each filter; the start and stop time and exposure 
for these observations are also given and (2) the best position 
derived from the GCN circulars. 

The Catalog contains a substantial amount of information 
related to the BAT instrument: (1) position, (2) significance of 
the detection, (3) type of trigger (image or rate), (4) length of the 
T90 and T50, (5) the start and stop of the T100, (6) T100 fluence 
in units of erg cnT 2 in the full energy range of 15-150 keV as 
well as four sub-ranges (15-25, 25-50, 50-100, 100-150 keV) 
and the spectral model used to estimate those quantities (power 
law or exponential cutoff power law), (7) 1 s peak flux in units 
of erg cnT 2 s -1 and of photons cnT 2 s -1 both in the full energy 
range of 15-150 keV and in the four sub-ranges, the time of the 
peak and the spectral model, (8) results of the spectral analysis of 
the T100 interval in the 15-150 keV band using the power law or 
the exponential cutoff power law model (spectral slope, energy 
of the cutoff, reduced y 2 , and degrees of freedom), and (9) 
hardness ratios of the fluence in the 25-50 keV and 50-100 keV 
ranges with respect to the fluence in the 15-25 keV. 


Those quantities are regularly published in the prelim- 
inary reports from the BAT team ground-based analysis 
(http://gcn.gsfc.nasa.gov/swift_gnd_ana.html). The Swift GRB 
Catalog simply collects them for an easier comparison and us- 
age. The BAT team, however, has run a more detailed and GRB- 
specific analysis. This work led to the publication of two catalogs 
(the First and the Second Swift BAT Gamma Ray Burst Catalogs, 
Sakamoto et al. 2008 and Sakamoto et al. 2011, respectively). 
The Swift GRB Catalog has been updated with the results from 
these peer reviewed analyses, in order to have more reliable and 
precise information. 

The Swift GRB Catalog also contains the following informa- 
tion from the GCN circulars or the literature: (1) the redshift 
(value, error, and type, i.e., if from absorption and/or from 
emission lines, or from photometric measurements), (2) the 
detection using ground-based telescopes (radio, infrared, and 
optical), (3) the trigger or follow-up observations with other 
X-ray and y-ray satellites, (4) the presence of a supernova sig- 
nature, (5) the association of the afterglow with a host galaxy 
(position, type of galaxy, redshift, and distance from the after- 
glow). 
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GRB09061 8: UVOT, TARG_ID=355083, Highest S/N Image GRB090618: UVOT, TARG_ID=355083, Finding Chart Image 

OBS_ID=00355083000, FILTER=WHITE OBSJD=00355083000, FILTER=V, EXPID=267006625 



Figure 11. Left Panel: UVOT FOV centered on the GRB using the summed image with the highest signal-to-noise ratio. The blue solid circle and the blue dashed 
annulus indicate the GRB and background-extraction region used to obtain the light curve analysis. The red dashed circle indicates the excluded field sources. Right 
panel: UVOT finding chart image with plotted BAT and XRT position error circles. Note the presence of trails for each detected source, indicating that Swift was still 
settling when those data were recorded. 

(A color version of this figure is available in the online journal.) 


GRB090618: UVOT, TARG_ID=355083 Magnitude Light Curves GRB090618: UVOT, TARGJD=355083 Integrated Flux Light Curves 



blue = Wl, magenta = M2, yellow = W2, orange = WHITE. As mentioned in the text, the user should ignore the first data point in the V filter if the satellite was still 
settling during the exposure (see also the previous figure). 

(A color version of this figure is available in the online journal.) 


8. COMPARISON OF PIPELINE RESULTS WITH 
PUBLISHED CATALOGS 

In order to test the reliability of the pipeline, a comparison 
with previously published results is mandatory. In particular, we 
focused our attention on works where the aim was to produce a 
catalog for one of the instruments on board Swift, i.e., Sakamoto 
et al. (2008) for the BAT catalog, Evans et al. (2009) for the 
XRT catalog, and Roming et al. (2009) for the UVOT GRB 
catalog. A consistency check was performed on a sample of 
27 GRBs that encompasses a wide range of behaviors in each 


instrument. For all the instruments we selected both long and 
short GRBs, as well as very bright and barely detected ones, and 
those also observed in the earliest and in the latest part of Swift’s 
mission. We tested the pipeline for some “extreme” cases: for 
the BAT, a few GRBs have a detected precursor while others 
were observed with a small partial coding fraction; for the XRT, 
we selected observations with high flickering between the WT 
and PC modes, with strong pile-up in both modes, and with 
the dead columns crossing the PSF on the side or perfectly on 
center; for the UVOT, we analyzed the emission from GRBs in 
crowded areas. 
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Figure 13. Combined BAT (black) and XRT (WT mode: blue and PC mode: red) light curves in units of erg cm 2 s 1 . The BAT light curve is converted into the 
0.3-10 keV range (left panel) and the 2-10 keV range (right panel). 

(A color version of this figure is available in the online journal.) 


Since the BAT and UVOT published catalogs extend up to 
2007 June 16, a comparison with the most recent observations 
was only performed for the XRT data. 

We chose as good testing parameters the fluence in the 
15-25, 25-50, 50-100, 100-150 keV ranges for the BAT, the 
optical/UV magnitudes for the UVOT, and the light curve slope 
decay for the XRT. 

Notes on the selection criteria that were used are as follows: 

1 . To estimate the BAT fluence, we considered the flux in each 
energy range obtained by fitting the spectrum in Xspec and 
we divided it by the exposure time. 

2. While Roming et al. (2009) performed the photometry on 
the individual sky images, regardless if a detection or an 
upper limit was found, our pipeline combines the individual 
sky images if consecutive upper limits are obtained in 
order to lower the detection/upper limit threshold. For this 
reason, only the optical/UV magnitudes in the individual 
sky images that have a detection in both Roming’s and our 
pipeline are considered. 

3. For the XRT comparison we used an approach similar to 
the one described in Evans et al. (2007). We downloaded 
the count rate light curve for all the GRBs from the online 
repository that was created by the University of Leicester. 
We performed the light curve analysis, finding the values of 
the time breaks where the light curve shows a change with 
respect to a constant (power law) decay. For each phase of 
constant decay, we estimated the slope of the power law 
decay. The same quantity was calculated in the analysis 
of the light curves produced by our pipeline, but with the 
assumption that the position of the time breaks is the same 
as the one previously estimated. In this way, a difference in 
the slope decay does not depend on the position of the time 
break. 

For all the testing parameters we evaluated the offset (in a) 
between the values obtained by our pipeline and from the other 
public catalogs. The offset is defined as (val p i pe — val cat )/a tot , 
where a tot is the propagated standard deviation and is defined as 
(CTp ipe +a c 2 at ) -1 / 2 . The offset distributions for the BAT fluence, the 
UVOT optical/UV magnitude, and the XRT light curve slope 
are shown in Figure 14. 

The plots show the reliability of our pipeline: the distribution 
of the a for the three quantities are Gaussian-like, with the 


peak close to zero. The standard deviations are 0.14, 0.35, and 

0.53 for the BAT, UVOT, and XRT distributions, respectively, 
meaning that 95% of the values produced by our pipeline are 
within 0.27a, 0.69a, and 1.04a from the values reported in the 
public catalogs. 

The discrepancies between the Swift GRB Catalog and other 
published works can be due to a large variety of factors. We 
name the most relevant as follows: 

1 . Calibration databases and software version releases (all), 

2. Data reduction filtering (all), 

3. Source and background-extraction region sizes (XRT and 
UVOT), 

4. Count rate correction methods (XRT and UVOT), and 

5. Background position (UVOT). 

9. SUMMARY AND CONCLUSIONS 

We have developed a semi-automatic pipeline that performs 
the data reduction and analysis of the three instruments on board 
Swift. The pipeline has been used to generate a uniform and 
complete analysis of all the GRBs observed by Swift and it will 
continue to be used as long as the mission is in operation. For 
the very recent GRBs, and for the ones that will be discovered 
or observed in the near future, the pipeline has used and will 
use the updated release of the CALDB files. The software tools 
that swifthlprod calls are also updated when bugs are fixed 
or new improvements are included. This will potentially affect 
the uniformity of the output products. For these reasons, we 
are planning to regenerate those products when major changes 
occur in the CALDB files or as released in the software tools. 

As mentioned in Section 1, there are many places where 
the user can retrieve data products for the GRBs observed by 
Swift. While our database is very much in sync with the products 
available at the University of Leicester Web site for the XRT data 
analysis, additional files are available for the BAT, especially 
regarding the spectral analysis: the PHA fits files during each 
single phase (pre-slew, slew, and post-slew, if applicable) can 
be downloaded, together with the respective response files, so 
that the user can perform a customized analysis. The main 
difference regarding the UVOT data products with respect to 
what the user can find in the Swift Burst Ground- Analysis Web 
site or at the “First Swift UVOT GRB Afterglow catalog” is 
that our light curves cover the full list of GRBs observed by 
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Figure 14. Distributions of the values obtained by our pipeline for the BAT fluence (left), UVOT magnitude (center), and XRT light curve slope (right), compared to 
the values reported in public catalogs. 


Swift and have been obtained by: (1) carefully selecting specibc 
background and source extraction regions, (2) using field stars 
in the estimate of the aperture-correction, and (3) regenerating 
the level 2 products using updated CALDB files. 

The results of the pipeline run were combined with data from 
the GCN circulars, the literature, and the BAT ground-based 
analysis to populate the Swift GRB Catalog, whose purpose is 
to generate a “summary” of the characteristics of those GRBs. 
We believe that this catalog is a useful tool for the public since it 
is constantly updated (typically within 1-2 business days after a 
GRB has been triggered or observed) and it collects information 
not specifically related to the Swift data analysis. As mentioned 
in Section 7, the catalog provides information regarding the 
measured redshift, the presence of a host galaxy, and if the GRB 
was observed in other wavelengths (radio, infrared, or optical) 
and by other high energy missions. This can help the user in 
selecting only those GRBs that meet specific criteria. 

Another important aspect that was taken into account in de- 
veloping the pipeline, is its versatility: indeed, it can be of 
help to perform an analysis of very different types of point- 
like sources, such as AGNs, pulsars, etc. This is particularly 
true when the user wants to analyze XRT and UVOT data (for 
an analysis of the BAT data for sources that are not GRBs, we 
suggest using the batsurvey tool developed by the BAT team). 
Changing the parameter settings for either swifthlprod, 
xrtgrblc, and/or uvotgrblc, one can easily optimize the use 
of those scripts to match specific needs. 

The script that handles the data reduction and analysis 
(swifthlprod), the results for each of the observed GRB 


(database), and the Swift GRB catalog are available for down- 
loading and querying at the HEASARC Web site, while the 
single scripts that perform the XRT and UVOT analysis 
(xrtgrblc and uvotgrblc) are already part of the HEAsoft 
release. 
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